arXiv:1505.02536vl [cs.IT] 11 May 2015 


Efficient Symbol Detection for the FSO IM/DD 
System with Automatic and Adaptive Threshold 
Adjustment: The Multi-level PAM Case 


Tianyu Song and Pooi-Yuen Kam 
Department of Electrical and Computer Engineering 
National University of Singapore, Singapore 117583 
Email: {song.tianyu, elekampy}@nus.edu.sg 


Abstract —To detect M -ary pulse amplitude modulation signals 
reliably in an FSO communication system, the receiver requires 
accurate knowledge about the instantaneous channel attenuation 
on the signal. We derive here an optimum, symbol-by-symboi 
receiver that jointly estimates the attenuation with the help 
of past detected data symbols and detects the data symbols 
accordingly. Few pilot symbols are required, resulting in high 
spectral efficiency. Detection can be performed with a very low 
complexity. From both theoretical analysis and simulation, we 
show that as the number of the detected data symbols used 
for estimating the channel attenuation increases, the bit error 
probability of our receiver approaches that of detection with 
perfect channel knowledge. 


I. Introduction 


In most current free space optical (ESO) communication 
systems, for reasons of simplicity, intensity modulation with 
direct detection (IM/DD) is used. To improve the spectral 
efficiency, we consider M-ary pulse amplitude modulation 
(MPAM) here by presuming that the receiver photodetector 
will not be saturated by the impinging optical signal level 
that has the highest energy. The transmitter side constellation 


is shown in Pig. 1(a) where I denotes the minimum signal 
intensity distance and the transmitter average power is defined 
as P = Y jl- Inside the receiver, as shown in Pig. 2 

there is an integrator that integrates the photo current for each 
symbol period Tg. For the fcth symbol interval {{k—l)Ts, kTg), 
the received electrical signal r{k) is obtained by sampling the 
integrator at time t = kTg. Since the photo current can be 
assumed to be constant during the integration time, r{k) can 
be expressed as |[Tj 

r(fc) = y/l/TsRhIm{k)Ts +n{k), (1) 


where R is the responsivity of the photo detector and h denotes 
the instantaneous channel gain. The transmitted data symbol 
m{k) takes on any value from set {0,1, ...,M — 1} with 
equal probability, and Gray mapping of bits onto the levels 
is assumed. In practice, multiplying the received signal by 
the normalising basis 0o(/) = i/l/Tg, which is shown in 
Fig. 0 is not necessary. We use it here because we want to 
simplify later performance analysis by normalizing the discrete 
additive white Gaussian noise (AWGN) term n{k) such that 
E[n(i)n(j)] = 5ijN[j/2 conditioned on that the continuous 
noise term n{t) is an additive white Gaussian random process 
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(a) Transmitter Constellation 
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(b) Receiver Constellation 
Fig. 1. Constellations 
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Fig. 2. Receiver diagram 

with mean zero and two-sided spectral density No/2. By 
defining A = ^l/TgRhlTs as the instantaneous receiver- 
side electrical-domain minimum signal distance, r(fc) can be 
modelled as © 

r{k) = Am{k) + n{k). (2) 


Correspondingly, the receiver side constellation is shown in 
Fig. 1(b) Since atmospheric turbulence and pointing errors 
cause fluctuations in the intensity of the received signal, i.e., 
h is time-varying, therefore, A is time-varying and is modelled 
as A = 2dh where 2d denotes the minimum signal distance 
when h = 1, i.e., 2d = ^Jl/TgRITg, which is related to the 
average transmit power P, i.e.. 


2d = s/ljTgRITg = 


2y/TsRP 
M-1 ' 


(3) 


For the on-off keying (OOK) system that is considered as 
the simplest MPAM {M = 2) case, the turbulence and pointing 





















error induced fading is the major issue that degrades system 
performance. Thus, the accurate channel state information 
(CSI), i.e., the accurate instantaneous value of h, is required 
for properly adjusting the decision threshold. In previous stud¬ 
ies, frequent insertions of pilot symbols Q or co-propagating 
reference light Q is used to estimate the CSI. In this paper, 
when it comes to pilot symbols, we specifically refer to 
pilot symbols used for channel estimation. In practice, pilot 
symbols are also required for other purposes such as timing 
synchronization. Since pilot symbols do not carry data, the use 
of them results in a system spectral efficiency reduction and 
additional energy overheads. Therefore, we hope to minimize 
the amount of pilot symbols. The co-propagating reference 
light has the same problem, and besides, since it has a different 
frequency with the signal-bearing light, the estimation result 
may not be accurate. For MPAM systems (M > 2), to detect 
signals reliably, the requirement of accurate CSI is more 
stringent. The difficulty of accurate CSI acquisition is the 
main reason that limits the use of MPAM in practical systems. 
However, though these two above-mentioned methods can 
be extended to MPAM systems, they are not so suitable in 
practice. The receiver design problem remains challenged by 
the fact that the CSI is hard to acquire at the receiver side if 
very few pilot symbols and no co-propagating reference light 
are used. 

Since the channel coherence length Lc, defined as the 
number of consecutive data symbol intervals over which the 
channel gain h can be considered to be constant, is a very 
large number (> 10^) for multi-Gbps systems ||^, the data 
and the CSI can be jointly detected and estimated. In 0, for 
the 2-level PAM system, i.e., the OOK system, based on the 
generalized likelihood ratio test (GLRT) principle, we derived 
a maximum likelihood (ML) sequence detection receiver (the 
GLRT-MLSD receiver) that jointly detects the data sequence 
and estimates the unknown channel gain. In Q, we have ex¬ 
tended this GLRT-MLSD receiver to multi-level PAM systems. 
However, as we discussed in Q, the search complexity of the 
GLRT-MLSD receiver increases with the modulation order M 
almost quadratically. Therefore, in section II of this paper, we 
derive a simpler decision-feedback (DFB), symbol-by-symbol 
receiver, whose implementation complexity is independent of 
the modulation order M. It uses the most recently detected 
data symbols to help estimate the instantaneous channel gain. 
As the number of the detected data symbols used to estimate 
the channel increases, the bit error probability (BEP) of this 
receiver approaches the Genie Bound, which is defined as the 
BEP of detection with perfect CSI (PCSI). 

H. The Decision-Eeedback Receiver 
A. Receiver Design 

We use rh{k) to denote the decision on symbol m{k) and 
m(fc, L) = [Th{k — L + 1), m{k — L-\-l), ..., m{k)] to denote 
the decision on the subsequence m(fc, L) = [m{k — L + 1), 
m{k — L+1), ..., m{k)]. The GLRT-MLSD receiver 

rh(fc,L)=arg max A(m(fc, L)), (4) 

m{/c,L) 


which was derived in 0 for OOK systems and generalised in 
0 for multi-level PAM systems, has been shown to be able to 
achieve the Genie Bound. In 0, A(m(A:, L)) is the decision 
metric and is given as 


A(m(fc, L)) 


(r(fc,L) • m(fc,L))2 
l|m(fc,L)P 


(5) 


where r{k, L) is the received signal subsequence with length 
L, i.e., r{k, L) = [r{k — L + 1), r{k — L + 1), ..., r(fc)]. This 
GLRT-MLSD receiver continuously and implicitly performs 
ML estimation of the by 0 Eq. (16)]. Here, for any hy¬ 
pothesized subsequence m{k,L), A = 2dh can be estimated 
similarly by performing 


A(m(fc, L)) 


r{k, L) ■ m(/c, L) 
||m(A:,L)P 


(6) 


Our GLRT-MLSD receiver has to evaluate the decision 
metrics of all possible subsequences and choose the one with 
largest metric value as the decision. In this section, based on 
0, we develop a new simpler decision-feedback symbol-by¬ 
symbol receiver. We first assume that at each time k, all the 
past decisions before time k have been completed, and we use 
[rh(fc — 1, L),m{k)] to denote the subsequence whose first L 
elements are rh[k — L), rh{k — L + 1), ..., m(/c — 1) and last 
element is m(k). Thus, at time k, the detection rule given in 
0 can be considered as choosing the m{k) that maximizes 
the decision metric value of subsequence [m(fc— 1, L), m{k)], 
i.e., 

m{k) = argmaxA([m(fc — 1, L),m{k)]). (7) 

m{k) 


Since [r(fc — l,L),r{k)] = r{k,L + 1) is independent of the 
detection result, after substituting 0 into 0, 0 is equivalent 
to 


m{k) 

= argmin (||[r(/c - l,L),r(fc)]f - A([rh(fc - l,L),m{k)\)) 

m{k) 

= arg min ||[r(A: — l,L),r{fc)] 

m{k) 

— A([rh(fc — 1,L), m{k)])[ni{k — 1, L), m(fc)]|p. 

( 8 ) 


As discussed in section 2, the channel coherence length Lc 
is on the order of 10"^, which is a very large number. The 
channel gain h can be safely considered unchanged from time 
point fc — 1 to fc. Thus, we can use the ML estimation on 
A = 2hd based on rh(/c — 1, L), denoted by 


A(rh(fc — 1, L)) = 


r{k — 1,L) ■ ih{k — 1, L) 


(9) 


|!m(A: - 1,L)||2 

to approximate A([rh(fc — l,L),m{k)]), with 

very high accuracy. This enables us to further 
simplify the decision rule 0 to ( [T0| . Since term 

r(fc — 1, L) — A{m.{k — 1, L))rh(fc — 1, L) is 

independent of m{k), we eliminate it and simplify the 
decision rule as 


rh{k) = arg min(r(fc) — A(rh(A: — 1, L))m{k))‘^. (11) 

m{k) 









m(fc) = arg min 

m{k) 


[r{k — 1, L),r{k)] — A{m{k — 1, L))[m.{k — 1, L), m{k)] 


2 


: arg mm 

m{k) 


v{k — 1, L) — Aiihik — 1, L))m(fc — 1, L) 


r{k) — A(m(fc — 1, L))m(k) 


( 10 ) 


In principle, to decide based on ( [TT] l, one has to compute the 
value of {r{k) — A(rh(fc — l,L))TO(fc))^ for each m{k) S 
{0,1, 2,M — 1} and then choose the m(k) corresponding 
to the minimum (r(k) — A{m.{k — l,L))m{kyf‘ value. Since 
m{k) may take any value from set {0,1, 2,M—1} and there 
are totally M entities in {0,1, 2 ,M — 1}, for one symbol 
detection, M evaluations of (r(fc) — A(rh(fc—1, L))m(fc))^are 
required. Thus, the search complexity still increases with M 
linearly. 

To reduce the complexity, we can further simplify the 
decision rule •HD as 


where ^ is the corresponding AWGN term. Substituting 
= A{M — 1) + into ([T3|) and simplifying, we have 


Ah — A ■ 


Y^i=l '^i,k 

Lm{M-iy 


(14) 


Obviously, Ak is a Gaussian random variable and its mean and 
variance are the tme value of A and — (M-iy ’ respectively. 
Apparently, we have 


,. No 

2L^(M- 1)2 


= 0 . 


(15) 


I 0 , r{k) < 0 

m{k) = < + |J ,elsewhere . (12) 

[ M-1 ,r(fc) > (M-l)ifc 

where Ak is obtained from (|^ and [ J is the floor function. 
In contrast to implementing implementing ( [T^ is much 
simpler. First, we need to compare r{k) to 0 and {M — l)Ak. 
If r{k) < 0 or r{k) > (M — l)Ak, the decision is made. If 
0 < r{k) < [M —l)Ak, we can substitute r(k) into + 
and obtain the decision. Totally, the computational complexity 
is much lower and is independent of the modulation order M. 

After implementing ([T^, a reverse Gray mapping is then 
performed to recovery information bits. 


B. The Channel Estimator 


1) A selective-store strategy: If all the elements of the 
detected subsequence th(fc— 1,L) are zeros, the denominator 
of A{ih{k — 1, L)) given in Q is zero, resulting in a channel 
estimation failure. We propose here a selective store strategy 
(SSS): at each time k, we only store the Lm most recent 
received signals that have been detected to carry symbol M—1. 
Since all the signals correspond to the same data symbol 
M—1, store of detection results is unnecessary, and thus, 
the estimate of A reduces to 


Ak 


M-1 


L„(M-1)’ 


(13) 


where is defined as the Ah most recent received signal 

at time k that is detected to carry data symbol M — 1. In this 
way, a zero-denominator of the channel estimator will never 
occur, leading to completely avoiding the channel estimation 
failure. 

2) Performance of the DFB receiver with the selective-store 
strategy: We assume all the decisions before time k are cor¬ 
rect. Even in practice, decisions cannot be 100% correct, with 
a very low error probability, we can assume so when analysing 
system performance. In the later section that shows simulation 
results, we do not assume 100% correct decision feedback. 
Thus, with the assumption, we have = A{M —l)-\-ni^k. 


From 0 Sec. 14.1], we see that for a random variable X, if 
E(Ar — a)^ = 0, where a is a constant, the random variable 
X equals a with probability 1. Thus, since Var(/i(in(/c — 
1, Lyf))) = E(/i(iii(fc — 1, Lu,)) — hy, we have 

lim Ak = A (16) 

Lm—^OO 

with probability 1. Since the estimation of A can approach the 
true value, we can say that the BEP of our DFB receiver with 
SSS approaches the Genie Bound as Lm goes to infinity. 

3) The generalised selective-store strategy: A more general 
SSS is given as follows; We selectively store the most 
recent received signals and their corresponding detection re¬ 
sults with the criterion fh > a, where a takes value from set 
{1,M — 1}. Clearly, a is an important parameter and the 
SSS introduced in the previous two paragraphs is an extreme 
case with a = M — 1. If 1, both detected data 

and corresponding signals are required to be stored, and the 
required memory length is 2Ljn', if a = M — 1, since only 
signals that are detected to carry M—1 are required to 
be stored, the memory length is Lm and the best memory 
efficiency is achieved. Also, with a fixed value of Lm, a higher 
value of a leads to better performance. Another aspect to 
consider is to keep the observation window length L^, which 
is on the order of LmM/{M — a), far smaller than the channel 
coherence length Lc- Thus, to achieve the best performance 
with a fixed Lm value, we can choose the highest value of a 
that makes Lml{f — a) < Lc/lO- 

III. Performance Results and Discussions 
A. The Genie Bound 

The BEP of symmetric PAM over the AWGN channel has 
been studied and given in 0 Eq.’s (9) and (10)] in terms of 
the average energy per bit Eb and the AWGN spectral density 
No- It should be noticed that Q assumed an electrical digital 
communication system and thus Eb is the electrical-domain 
average energy per bit. Based on the mapping given in 0 
Eq. (2)] which is from Eb to the minimum signal distance 
2d, we can obtain its inverse. After substituting this inverse 
into 0 Eq.’s (9) and (10)], we can derive the general BEP 


















(17) 


jPAM 


M! 

Nn 


log2 M 

.")= i: 

k=l 




i=0 


M log2 M 


expression in terms of the minimum signal distance A and the 
AWGN one sided spectral density Nq as ([nil. Since © is 
expressed in terms A and Nq, it applies to both symmetric and 
asymmetric PAM signals over the AWGN channel. Therefore, 
the BEP of the PCSI receiver conditioned on a given value of 
h is ,M), and the average BEP over all possible 

values of h is given by 


n''“'(e)= r P!^'‘{^A^.M)n{h)dh. ( 18 ) 

Jo .'''0 

This average BEP, which is also referred to as the Genie 
Bound, is used as a benchmark when analysing other receivers. 
Specifically, when discussing the BEP of our receiver, we refer 
to the average BEP over all possible channel states. 

Prom ( frsj l, we see that the average BEP is only related to 
2d, Nq, M and Ph{h). Thus, when the environment condition, 
i.e., Ph(h), the system modulation order M and the receiver 
circuit AWGN one-sided PSD Nq do not change, using 
different transmission rates (different bandwidths), to achieve 
the same error probability, we need to keep, d unchanged. 
Por conventional electrical digital communication systems, we 
have 


M-l 


— i^{2df = Eblog^M, 

i=0 

which is equivalent to 


2d = 


' 6Eb log2 M 
(M- 1)(2M- 1)’ 


(19) 


( 20 ) 


where E^ denotes the electrical-domain energy per bit. From 
we see that when keeping the same level of error proba¬ 
bility, Eb is irrelevant to the transmission rate. This means no 
matter how fast we transmit data, the energy consumed on each 
bit does not change, i.e., to transmit the same amount of data, 
the total energy consumed is irrelevant to the transmission rate. 

However, for optical systems, the situation is different. The 
optical-domain energy consumed per bit E^ can be calculated 
via 


E°b = 


PTs 

log2 M ■ 


According to ([^, we have 


^2 ( V%Rp V {Ej) 


\M -I ) 


T. 




P? log 2 M 
(M- 1)2 ’ 


Eiog^y 

M-l ) 


( 21 ) 


( 22 ) 


where i?data = (log 2 M)/Ts is defined as the data rate. We can 
see that the higher the transmission rate i?data is, the lower the 
energy consumed on each bit transmission is. Therefore, we 
can say that, when 2d, Nq, M and ph{h) are unchanged, i.e., 
the same level of error probability is retained, to transmit the 
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i ■ 2'=-! 
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Q 



same amount of data, using a larger data rate (bandwidth) is 
suggested since the total energy consumption is lower. 

This is because the amplitude of the photo detector output 
current is proportional to the incident light power Q. That 
means if we increase the optical power by S times, the 
electrical-domain signal amplitude is increased by S times and 
the electrical-domain signal power is increased by times. 
When we use a larger bandwidth, say K times of before, 
to achieve a higher data rate, the noise power, which is the 
product of the noise power spectral density and the bandwidth, 
is increased by K times. To keep the error probability same, 
the electrical-domain signal power is also required to be 
increased by K times, corresponding to increasing the optical 
signal power by y/K times. Since the transmission rate is 
increased by K times, the required time for transmitting the 
same amount of data is shortened to ^ of before, and thus the 
total optical energy consumption is ^/K x P = times of 
before. Finally, our conclusion is that using a higher bandwidth 
is suggested since it saves energy. 

It should be emphasized that all the analysis and conclusions 
in this subsection are based on several basic assumptions; 


1. pin diode is used for photo detection, i.e., thermal noise 
can be regarded as the dominant noise souse and shot 
noise is negligible 

2. the photo detector will never be saturated by the signal 
pulse that is with the highest power level 

3. the photo detector can response as fast as the user 
requires, i.e., the photo detector bandwidth is no less than 
the system symbol rate. 

4. the FSO channel is always in a linear regime, where the 
output power is proportional to the input signal power. 


Without any of the abovementioned assumptions, the analysis 
and conclusions may not be true. 


B. Simulation Parameters 

As discussed in Q and eg, geometric spread and pointing 
errors hp, atmospheric turbulence ha, and path loss hi together 
determine the overall channel gain h. The channel state h 
is formulated as h = hphahi. All these three parameters 
are time-varying for outdoor environments. The path loss is 
related to the atmospheric phenomena, such as rain, fog and 
snow, and the distance between the transmitter and receiver. 
Therefore, if the distance between the transmitter and the 
receiver does not chain, hi changes much slower than ha 
and hp. In Q and | [T0| , hi is considered as a deterministic 
variable while ha and hp are considered as random variable. 
In fTT] , log-normal distribution is adopted to model ha for 
weak turbulence, Gamma-Gamma distribution for moderate to 
strong turbulence and the negative exponential distribution for 
strong turbulence. Since in | [T2) , it has been shown that the 
Gamma-Gamma distribution can nicely fit the channel fading 
statistics of all turbulence regimes, in this paper, we only 


















consider ha is a Gamma-Gamma distributed random variable, 
and the pdf of ha is 




2(q,^)(«+/3)/2 

r(a)r(/3) 


^(a+/3)/2-l 


-^CK—/5 


(2^/aPh^ 


h > 0, 
(23) 


where Ka{-) is the modified Bessel function of the second 
kind, and 1//3 and 1 /a are the variances of the small and 
large scale eddies, respectively. Pointing error influence on an 
FSO system is discussed in | |T0| and |13|, and we here use 
the model in pO) where the pdf of hp is given as 


Php{h) = \ 0<h< Aq. 


(24) 


Parameter Aq is the fraction of the collected power when 
no pointing etTor occurs, and 7 is the ratio between the 
equivalent beam radius at the receiver and the pointing error 
displacement standard deviation at the receiver p0| . Without 
loss of generality, we can incorporate hi into ha which 
amounts to setting hi = 1. Then for a turbulent channel with 
pointing errors, the channel gain is h = hahp, and by using 
||^ Example 13.1.10], its pdf can be derived by 

Ph{h) = J da, h > 0. (25) 

Full model details can be found in 0 , 0 , (D-jig. 

The system parameters of this paper are given as follows: for 
the weak turbulence channel, a = 17.13, /3 = 16.04 and the 
corresponding SI = 0.1244; for the strong turbulence channel, 
a = 2.23, /3 = 1.54 and SI = 1.3890. According to |[T4), 
the turbulence with SI=0.1244, which is less than 1, is in 
weak irradiance fluctuations regime; and the turbulence with 
SI= 1.3890, which is larger than 1, is in moderate-to-strong 
irradiance fluctuations regime. The pointing error parameters 
are chosen to be Aq = 0.0198 and 7 = 2.8071. Without loss of 
generality, the photo detector responsivity R is assumed to be 1 
and the expected channel gain E[/i] is set to be 1. Additionally, 
we assume the system data rate i?data is on the order of 
10 Gbps, thus the symbol duration Tg = (log 2 M)/i?(iata is 
on the order of 10 “® — 10 “^°s and the channel coherence 
length La can be safely regarded as 10^ data symbols. We 
consider the typical thermal noise, which is -174dBm/Hz, 
passing through a 5012 receiver circuit. Thus, the value of 
No/2 is -174dBm/Hz ^ 5012 = 10“2O-4w/Hz -7 5012 = 
7.96 X 10“23AVHz, i.e., Nq = 1.59 x lO-^^A^/Hz. 


C. Numerical Results and Discussion 

In Fig. 0 we plot the BEP versus average receive power 
curves of our DFB receiver. Since the implementation com¬ 
plexity of our DFB receiver is independent of M, we show 
simulation results with up to M = 32. We see that as 
the value of increases, the BEP decreases. To achieve 
the Genie Bound, the required value of is 12 for the 
weak turbulence channel and no larger than 16 for the strong 
turbulence channel. 



Average Receive Power (dBm) 
(a) SI = 0.1244. 



Average Receive Power (dBm) 

(b) SI=1.3890. 

Fig. 3. Performance of M-PAM signals. 

From Fig. we see that if we keep the data rate unchanged 
and increase M to save bandwidth, to achieve the same error 
probability, we need to increase the transmit power. This 
means bandwidth and power are a pair of contradictions. We 
can hardly minimize both of them simultaneously, but we can 
select appropriate values according to practical requirements. 

Another interesting observation is that if we keep M un¬ 
changed and use a larger bandwidth to achieve a higher data 
rate, to achieve the same level error probability, the multiple 
of the power growth is the square root of that of the data 
rate growth. Specifically, from Fig. 0 we can see that for 
M = 4, if we increase the data rate from 10 Gbps to 20 
Gbps, the corresponding power increment is less than 3dB, 
approximately 1.5 dB. 

To further study the receiver performance with more values 
of memory lengths, we plot Fig. 0 For the reason of space 
limit, we only pick one {M, average receive power, data 
rate, SI) point here, which is (M = 16, -1 dBm, 40 Gbps, 
SI=1.3890). Clearly, as the value of Lm increases, all the BEP 











IV. Conclusions 



Since the channel coherence length Lc is very large, we can 
use the detected data symbols to estimate the unknown channel 
state instead of using pilot symbols which causes spectral 
efficiency reduction. In this paper, based on the decision 
metric of our previously proposed GLRT-MLSD receiver, we 
propose a DFB symbol-by-symbol receiver, whose implemen¬ 
tation complexity is much lower and is independent of the 
modulation order. Hence, it is efficient both spectrally and 
computationally. We also propose a selective-store strategy, 
which can help avoid potential channel estimation failures and 
increase the system memory efficiency. Additionally, we derive 
a general BEP expression for both symmetric and asymmetric 
M-PAM signals, by using which as a benchmark, we have 
shown that as the number of the detected data symbols used 
to estimate the channel increases, the BEP of our DEB receiver 
approaches the Genie Bound. 


Fig. 4. Performance with different Lm’s; M = 16, -1 dBm, 40 Gbps, 
SI=1.3890. 
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(a) Lm — 1. 




(c) Lm = 8. 


Fig. 5. Estimates of I with different Lm’s', M = 2, ^ = 20dB 


curve approach the corresponding Genie Bound. 

In Eig. 1^ we plot the estimate of A with different memory 
lengths. We assume OOK modulation and choose ^ = 20dB. 
The numerical results given in Fig. completely agree with 
the theoretical analysis given in the previous section. Thus, we 
have shown that, by both theoretical analysis and simulation, 
with a higher value of Lm, the estimation of A is more accu¬ 
rate. This enables the BEP of our DEB receiver to approach 
the Genie Bound. 

When starting to operate our DEB receiver, Lm pilot sym¬ 
bols are required to obtain an initial value of I. Unlike some 
other communication systems that require frequent insertion 
of pilot symbols to avoid burst errors, in our simulation, no 
burst errors are observed even with only Lm pilot symbols for 
initializing and no further frequent insertions. Thus, the pilot- 
to-data ratio is far smaller than 1, and we do not consider the 
power consumed by pilots when calculating the average SNR 
per bit. 
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